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Abstract 

We obtain an asymptotic expansion for the sum 


S{a-w) = ^ 

71=1 


e~ an 


2 


n w 


as a —y 0 in | arg a\ < for arbitrary finite w > 0. The result when w = 2m, where m is 
a positive integer, is the analogue of the well-known Poisson-Jacobi transformation for the 
sum with m = 0. Numerical results are given to illustrate the accuracy of the expansion. 
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1. Introduction 


The classical Poisson-Jacobi transformation is given by 


n= 1 



(i.i) 


where the parameter a satisfies 5R(a) > 0. This transformation relates a sum of Gaussian 
exponentials involving the parameter a to a similar sum with parameter ir 2 /a. In the case 
a -A 0 in Jf(a) > 0, the convergence of the sum on the left-hand side becomes slow, whereas the 
sum on the right-hand side converges rapidly in this limit. Various proofs of the well-known 
result (1.1) exist in the literature; see, for example, [3, p. 120], [4, p. 60] and [5, p. 124]. 

In this note we consider the sum 


OO 

%;») = X 

n= 1 


n w 


(5ft(a) > 0). 


( 1 . 2 ) 


This sum converges for any finite value of the parameter w provided 5R(a) > 0; when a = 0 
then S(0;w) reduces to the Riemann zeta function C(w) when Jf(u;) > 1. Consequently, the 
series in (1.2) can be viewed as a smoothed Dirichlet series for C( w )- The asymptotic expansion 
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of S(a;w) as a —> 0 in 5ft(a) > 0 is straightforward. The most interesting case arises when 
w = 2m, where m is a positive integer, for which we establish a transformation for S(a;2m) 
analogous to that in (1.1) valid as a —> 0 in 5ft(a) > 0. This similarly involves the series in 
(1.2) with a replaced by 7r 2 /a, but with each term decorated by an asymptotic series in a. A 
recent application of the series with w = 2 and w = 4 has arisen in the geological problem of 
thermochronometry in spherical geometry [6]. 


2. An expansion for S(a; w) as a —> 0 when w ^ 2,4,... 

Our starting point is the well-known Cahen-Mellin integral (see, for example, [3, §3.3.1]) 


2 /*C+OOi 

= — / T(s - a)z~ s ds (z ^0, | arg z\ < \tv), 

Z7TZ Jq—qqI 


( 2 . 1 ) 


where c > 5ft(a) so that the integration path passes to the right of all the poles of T(s — a) 
situated at s + a — k (k = 0,1, 2,...). For simplicity in presentation we shall assume throughout 
real values of w > 0. Then, it follows that 


p —an ^ n~ w r c +°° l 

S(a-w) = ]T — ~ = / r(s)(an 2 ) _s ds 

n=1 n—1 J c-ooi 

rC+OOi 

— / r(s)C{2s + w)a~ s ds, 

J c—ooi 
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upon reversal of the order of summation and integration, which is justified when c > max{0, ^ — 
and evaluation of the inner sum in terms of the Riemann zeta function. The integrand 
possesses simple poles at s = 1 and s — — k ( k = 0,1,2,...), except if w = 2m + 1 is an 

odd positive integer when the pole at s = ^ is double. The case when w = 2m is an even 

positive integer requires a separate investigation which is discussed in Section 3. 

Consider the integral taken round the rectangular contour with vertices at c±*T, —c'HT, 
where c 1 > 0. The contribution from the upper and lower sides s = a ± iT, —d < a < c, 
vanishes as T —> oo provided | arg a| < since from the behaviour 

F(a± it) = 0(t a ~h-^ t ), C(v±it) = 0(t^log A t), {t oo), 


where for cr and t real 


fi(a) = 0 (ct > 1), \ — ^cr (0 < a < 1), \ — a (a < 0), 

A = 1 (0 < a < 1), A = 0 otherwise, 

the modulus of the integrand is controlled by 0 (T 0 ' + A <T ) _ 5 log Te _AT ), with A = 5 7r — I arg a|. 
The residue at the double pole s = —m when w = 2m + 1 (m = 0,1,2,...) is given by 

( _ n ) m 

-j —{7 - \ log a + \i)(m + 1 )}, 

m! 

where 7 is Euler’s constant and ip(x) is the logarithmic derivative of the gamma function. 
Displacement of the integration path to the left over the poles then yields (provided w ^ 2m) 

N—1 / 

S(a ; w) = J(a; w) + ((w — 2 k)a k + Rn , 

/ Id 
k =0 


( 2 . 2 ) 
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where 


i r a - i w ) a ^- L )/ 2 


(w ^ 2TO + 1) 


J(a; w) = 


(-ay 


ml 


'{7 - \ log a + hy , ( m + 1 )} (w = 2 ?n+l), 


iV is a positive integer such that N > + i and the prime on the sum over k denotes the 

omission of the term corresponding to k = m when w = 2to + 1 . 

The remainder Rjy is 


-1 /»—c+ooi 

Rn = -— t / r(s)£(«; + 2s)a~ s ds, 

ZiTTl J _ c — QQ2 


c = N - 


(2.3) 


It is shown in the appendix, when w y 2,4,..., that = 0(a N ~ 3 ) as a —>• 0 in | arg a\ < |7r, 
with the constant implied in the O-symbol growing at least like T(lV-)-l—^tu). This establishes 
that the above series over k diverges as N —» 00 and that (2.2) is therefore an asymptotic 
expansion. 

We remark that the algebraic expansion (2.2) also contains a subdominant exponentially 
small component as a —> 0; compare [3, §8.1.5] for the particular case w = 0. We do not 
consider this further in the present paper. 


3. An expansion for S(a\ 2m) when m = 1,2,... 


The case w = 2m, where m is a positive integer, is more interesting as this leads to the 
analogue of the Poisson-Jacobi transformation (1.1). There is now only a finite set of poles of 
the integrand in (2.1) at s = \ — |w and s = 0, —1, —2,..., — m, since the poles of T(s) at 
s = —m — k (k = 1,2,...) are cancelled by the trivial zeros of the zeta function ((2m + 2s) 
at s = —m — 1, —m — 2,.... This has the consequence that the integrand is holomorphic in 
Sft(s) < —m, so that further displacement of the contour can produce no additional algebraic 
terms in the expansion of S(a\ 2m). Thus, we find when w = 2m 


m f \k 

S(a ; 2to) = |r(| - m)a m ~2 + ^ — ((2m — 2k) a k + II, (3.1) 

k-0 


where, upon making the change of variable s —> — s, 

II = - [ T(—s)((2m — 2s)a s ds 

2 t n J L 

and L denotes a path parallel to the imaginary axis with Jf(s) > to. 
We now employ the functional relation for ((s) given by [5, p. 269] 

((s) = 2 s 7t s-1 £(1 — s)T(l — s) sin hns 


(3.2) 


(3.3) 


to convert the argument of the zeta function in (3.2) into one with real part greater than unity. 
The integral in (3.2) can then be written in the form 


(-) m (27r) 2m 

2-7U 



2 TO + 1) 


r(2s — 2 to + 1) 

r(s + i) 



ds. 


Since on the integration path 5R(2s — 2 to+ 1) > 1, we can expand the zeta function and reverse 
the order of summation and integration to obtain 


II = R m 7r 2m -5 Y, n 2m " 1 K n (a-,m), (3.4) 

n —1 
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where 


K n (a; to) := 


i f r(s — to + i)r(s — to + 1 ) 


(M 


2m J L r(s + 1) 

and we have employed the duplication formula for the gamma function 

r(2 z ) = 2 2 z_ 1 7 t“^ r(*)r(* + \). 


ds. 


The integrals K n (a: to) have no poles in the half-plane 5ft(s) > to, so that we can displace 
the path L as far to the right as we please. On such a displaced path |s| is everywhere large. 
The quotient of gamma functions may then be expanded by making use of the result given in 
[3, P- 53] 


r(* m+ ^ T ^ s w + 1) = ^ (- y Cj T(a + 0-j) + Pm{s)T(s + ■& - M) (3.5) 


r(s +1) 


3=0 


for positive integer M, where d = \ — 2m, 

_ (w )j{m+\)j _ 2~ 2j '(2to) 2 j 

C ° j! f- 

and pm{s) = 0(1) as |s| —> oo in | arg s| < n. Substitution of this expansion into the integrals 
K n (a;m) then produces 

M — l 


K n (a; to) = + (^ 2 ) ds + K M 


3=0 

M—l 


I y j ( a y 3 2 -n 2 n 2 /a _, -p 

2s (~) C 3 e +1Zm 


(3.6) 


3=0 


by (2.1), where 


TZm = E J pm(s)T(s + d- M) {-^ 2 ) ds - 


Bounds for the remainder 7 Zm have been considered in [3, p. 71, Lemma 2.7], where it is shown 
that 


Um °(( 7 r 2 n 2 ) 


a \M— 1 ) _ 22 


Z Z / 

—tv n /a 


(3.7) 


as a —> 0 in the sector | arg a| < 

Collecting together the results in (3.2), (3.4), (3.6) and (3.7), we obtain 




1 OO __ 7r 2 r? 2y 0 | M — l 


-,2m 


E 

3=0 


—a 


O 


(s?) 


M 


From (3.1) we now have the following theorem: 

Theorem 1. Let to and M be positive integers. Then, when w = 2 to, we have the expansion 
valid as a —> 0 in | arg a\ < 


m / \ ^ 

S(a; 2 to) = |r(^ — m)a m ~i ^ C(2m — 2k) a k 

k—0 
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where T n (a;m) has the asymptotic expansion 


T n (g; to) __^ n 2 /o 


T„(a;m) = 2^ -^- 

J=o J ' 


(4*V 

\7r z n z / 


This is the analogue of the Poisson-Jacobi transformation in (1.1). In the case m = 0, the 
quotient of gamma functions in (3.5) is replaced by the single gamma function T(s + 1), with 
the result that cq = 1, Cj = 0 (j > 1) and T„(a;m) = 1 for all n > 1. Then (3.8) reduces to 
( 1 . 1 ) and is valid for all values of the parameter a (not just a —> 0 ) satisfying | arg a| < 

Remark 1. We note that the values of the zeta function appearing in (3.8) can be expressed 
alternatively in terms of Bernoulli numbers by the result [2, p. 605] 

(2ir) 2n 


4. Numerical results and concluding remarks 

From the well-known values [2, p. 605] 

we obtain from Theorem 1 the expansions in the cases m = 1 and m = 2 given by 


_ , n a . .i /a \2 




0(a M ) 


A 2 

CU A\ V n ~ a 

S(a: 4) =- 

v ; 90 6 


2 1 3 

— 7 r 2 a 2 

3 


7 OO ___2_2 / 


O’i: 


f W (!)j(2)j 


0(a M ) 


valid as a —>■ 0 in | arg a] < ^tt. 

In Table 1 we show the results of numerical calculations for the case m = 2. For different 
values of the parameter a we present the value of the absolute error in the computation of 
iS(a;4) from (4.2). In the computations, we have used only the n = 1 term (since the order of 
2 exp(— 47 r 2 /a) was found to be less than the error), with the expansion for Ti(a; 2 ) optimally 
truncated (corresponding to truncation at, or near, the least term in modulus) at index jo ~ 
( 7 r 2 /a) — |. It is seen that the error when a = 0.1 is extremely small and that, only when a ~ 2 
does the relative error start to become significant. 

To conclude, we mention that a similar treatment can be carried out for the sum 


S p (a; w] 


= £ 


(a -»■ 0 , S(a) > 0 ) 
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Table 1: Values of the absolute error in the computation of ,S'(a; 4) from (4.2). The value of the index jo 
corresponds to optimal truncation of the expansion Ti(a; 2). 


a 

5(o;4) 

Error 

Jo 

0.10 

0.952696 

9.662 x 10~ 86 

96 

0.20 

0.849025 

9.768 x 10~ 43 

46 

0.25 

0.803169 

4.045 x 10~ 34 

36 

0.50 

0.615128 

7.769 x 10- 17 

17 

0.75 

0.475493 

4.656 x 10- 11 

10 

1.00 

0.369026 

3.642 x 10~ 8 

6 

1.50 

0.223285 

2.856 x 10~ 5 

3 

2.00 

0.135356 

7.500 x 10- 4 

1 


for positive even integer w and p. The case w = 0 and p > 0, corresponding to the Euler-Jacobi 
series, has been considered in [3, §8.1]; see also [1] for a hypergeometric approach when p is a 
rational fraction. The details of the small-a expansion of S p (a;w) will be presented elsewhere. 


Acknowledgement. The author wishes to acknowledge B. Guralnik for having brought this 
problem to his attention 


Appendix: A bound for the remainder Rjy 


Let ip = arg a and integer N > \vj + Upon replacement of s by —s followed by use of (3.3), 
the remainder Rn in (2.3) becomes 


Rn = 


(27 


nN— i+002 


C(1 -w + 2s) 


2m J N _i _ooi ' T(1 + s) 

With s = N — i + it, t £ (—oo, oo) we have 


T(1 - w + 2s) sin7r(s - \w) ( a y 


ds. 


sin tts 


l«»l s I 2 ')” -1 (i) 


N-i 


C(2 N - w) 


<!>t 


Y[2N - w + 2 it) 


T{N + ^ + it) 


dt , 


since |()(:r + if)| < ((x) (x > 1) and 


sin7r(A — \ + it) 


sin7r(A + \ + it) 

It then follows that 

jV-i roc 


| 008 77(111; — it) | (cos 2 4Tin + sinh nt) 2 


cosh 7 rf 


\R N \=o((±) e~^\T(N-lw + it)\ 


cosh 7 it 


T (N — + 1 + it) 


< 1. 


T(N + \ + it) 


dt . 


(A.l) 


Using the argument presented in [3, p. 126], we set N — 1 w — \ = M + 5, with — \ < <5 < \ 
so that M < N — 1, to find 


T(N — \ + it) 


T(N- 


it) 


= P(t)g(t), g(t) := 


r(i + <5 + it) 


T(1 + it) 
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where 


pit) = a + o-i nJhFi+V+eii 

From the upper bound for the gamma function r(z) with 2 
|r(2)| < T(x){l + t 2 /x 2 )^ x -ie-^ Wt) e 1/{6 ' zl \ 


<(± + t 2 )-i <2. 

= x + it, x > 0 [3, p. 35] 
cj)(t) = arctan(|t|/x) 


< r(s)(l + t 2 )^-3 exp [*{w(r) - ±tt|t|} e 1/(6a:) 

< e x T(x)(l + T 2 )2*-J e -3™M e 1 /^, 

where we have put r = t/x, defined w(r) = |r| arctan(l/|r|) and used the fact that 0 < w(r) < 1 
for r € [0,oo), with the limit 1 being approached as t —> 00. Substituting the above bounds 
into (A.l), we see on setting x = N — that 

\R n \ = e w r(7V-±u> + l)0 2 J (1 + T 2 ) Ar/2 g(r){e _A+r + e _A ~ r } dr ^ , 

where A± = (N — iw)(i7r ± ^>). Since g(r) = 0(r <5+ ^) as r —>■ 00, the integral is convergent 
provided |t/;| < ^71 and is manifestly an increasing function of N. 

Hence 

Rn = 0(ci n ~^) (a0, |arga|<i7r), (A.2) 

with the constant implied in the O-symbol growing at least like r(N + l — ^w) as N increases. 
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